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Channeling and blocking in crystals refer to the orientation dependence of ion penetration in 
crystals. "Channeling" occurs when ions incident upon a crystal along symmetry axes and planes 
suffer a series of small-angle scatterings that maintain them in the open "channels" in between rows 
or planes of lattice atoms and thus penetrate much further than non-channeled ions. Channeled ions 
do not get close to lattice sites, where they would suffer a large-angle scattering which would take 
them out of the channel. "Blocking" is the reduction of the flux of ions originating in lattice sites 
along symmetry axes and planes due to the shadowing effect of the lattice atoms directly in front 
of the emitting lattice site (see e.g. the review by D. Gemmell and references therein). 

Channeling and blocking effects in crystals are extensively used in crystallography, in mea- 
surements of short nuclear lifetimes, in the production of polarized beams etc. Channeling must 
be avoided in the implantation of B, P and As atoms in Si crystals to make circuits. The channel- 
ing effect in Nal(Tl) crystals was first observed by Altman et al. in 1973 who observed that 
channelled ions produce more scintillation light because they loose practically all their energy via 
electronic stopping rather than nuclear stopping. When ions recoiling after a collision with a WIMP 
move along crystal axes and planes, they give all their energy to electrons, so the quenching factor 
Q (the fraction of the energy deposited in a collision that goes into scintillation or ionization) is 
2 ~ 1 instead of e.g. Qi ~ 0.09 and Qj^a — 0.3 in Nal(Tl). The potential importance of this effect 
for direct dark matter detection was first pointed out for Nal(Tl) by Drobyshevski in 2007 and soon 
after by the DAMA collaboration [^. This last paper gave an estimate of the channeling fraction 
of recoiling ions as function of the energy, in which the fraction grew with decreasing energy to 
be ~ 1 close to 1 keV. With this channeling fraction, the region of light WIMPs compatible with 
producing an annual modulation as measured by the DAMA collaboration shifted considerably to 
lower cross sections, by about one order of magnitude (see e.g. [Q]). Channeling could also produce 
a novel dark matter signature. Since the WIMP wind comes preferentially from the direction fixed 
to the galaxy towards which the Sun moves, Earth's daily rotation makes this direction change with 
respect to the crystal, which could produce a daily modulation in the measured energy (equivalent 
to a modulation of the quenching factor). As first pointed out by Avignone, Creswick, Nussinov in 
2008 [||], this modulation depends on the orientation of the crystal with respect to the galaxy, thus 
it would be a background free dark matter signature. 

Our calculations [^] of the fraction of recoils that are channeled as function of recoil energy 
and direction use classical analytic models which started to be developed in the 60's and 70's, 
soon after channeling was discovered, for ions of energy MeV and higher. We use in particular 
Lindhard's model [^, supplemented by the planar channel models of Morgan and Van Vliet 
and the 1995-1996 work of G. Hobler on low energy channeling (applied to ion implantation 
in Si). In these models the rows and planes of lattice atoms are replaced by continuum strings 
and planes in which the screened Thomas-Fermi potential is averaged over a direction parallel to a 
row or plane of lattice atoms to find the continuous potential U (see Fig. l.a for examples of U). 
The appearance of continuous strings or planes can be understood as the overlap of the "Coulomb 
shadows" of individual atoms in a lattice row or plane behind the direction of arrival of a parallel 
beam of positive ions, when the beam arrives at a small enough angle, smaller than a critical angle 

Then, the individual shadows overlap forming a string or plane of thickness Pc within which 
the scattered ions do not penetrate (see e.g. Fig. 4.6 of Ref. [p^). We included just one string or 
plane, which is a good approximation except when Pc approaches the radius of the channel. 
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Figure 1: a. (left) Continuum potential for a string (black) and plane (green) for a Si ion propagating in Si 
as function of the distance perpendicular to the string or plane, b. (right) Oberved (green) and predicted with 
c = 2 (black) channehng angles for B ions propagating in Si at room T. 

Lindhard proved that for channeled ions the "transverse energy" Sperp = + ?7 ~ 

E^'^ + ?7 is conserved, where M, E = Mv^/2 are the mass and kinetic energy of the propagat- 
ing ion, Vperp = vsin(/> ~ v(^ is the component of the velocity perpendicular to the string or plane 
and U is the continuum axial or planar potential at the position of the ion. Thus, for a particu- 
lar fiperp (determined by the recoil energy E, initial recoil angle 0,- and initial potential and Ui), 
a channeled propagating ion does not approach the string or plane closer than a minimum dis- 
tance Pmin (for which Vperp = 0) and far away from the string or plane, close to the middle of 
the channel, the ion moves on a trajectory forming an angle ^|f with the string or plane, given by 
Sperp = E(^f + Ui = U (pmin) = E^^ + ^^middie- Channeling occurs when py^m > Pc> which amounts 

tOY= \/[t/(Pmin)- ^middle)] < ¥c = V P (Pc) - ^ middle)] / E . 

All the difficulty of this approach resides in calculating Pc. Including temperature (T) effects. 



the critical distance can be approximated by Pc{E, T) = y Pc{E) + [c ui (r)]^, where Pc{E) is the 
critical distance for a perfectly rigid lattice (which decreases as E increases, as shown in Fig. 2. a), 
ui{T) is the 1 -dimensional amplitude of thermal fluctuations (we used the Debye model) which 
increases with T, and c is a dimensionless number which (for several crystals and propagating 
ions, different than the ones we study) was found through data and simulations to be 1 < c < 2. At 
large enough energies Pc{E,T) ~ c ui{T) and thus as T increases the strings are planes become 
thicker, the channels narrower and y/^. smaller (see Fig. 2). Using this formalism with c = 2 we 
could reproduce data on channeling angles of B and P ions in Si measured at room T provided by 
Hobler |^ (shown in green for B in Fig.l.b). As shown in Figs, l.b and 2.b the critical angle Yc 
increases with decreasing energy until the critical distance for channeling p^ approaches the radius 
of the channel, at which point \j/c goes zero (see Fig. l.a). At still lower energies, Pc should be 
larger than the radius of the channel, what indicates that channeling is not possible (thus Xj/c = 0). 

The channeling of an ion depends not only on the initial angle its trajectory makes with a 
string or plane, but also on its initial position. The DAMA collaboration calculated the chan- 
neling fraction for ions starting their motion close to the middle of the channel, where channeling 
happens if y < ¥c- However this is not the case in dark matter detectors, since the recoiling 
ions start their motion at or close to their original lattice sites (and leave those sites empty). For 
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Figure 2: a. (left) Critical channeling distance and b. (right) critical channeling angles \j/c for a Si ion 
propagating in the <100> channel of a Si crystal, for a static lattice (dashed green) and T corrected (solid 
black lines) at different T with c = 2 (2m i shown in red and the radius of the channel in solid green) 
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Figure 3: Channeling fractions for a Si ions ejected from lattice sites in a Si crystal as function of the ion 
energy for different temperatures T and T effects computed with a. (left) c = 1 and b. (right) c = 2. 
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Figure 4: Same as Fig. 3 but for Ge ions ejected from lattice sites in a Ge crystal. 
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Figure 5: a. (left) Upper bounds to the channeling fraction of Na and I ions ejected from lattice sites in a Nal 
crystal for several temperatures (with c — 1). b. (right) Using the T= 293 K curve as the channeUng fraction 
in the DAMA detectors the region of compatibility with the DAMA signal for WIMPs with spin-independent 
interaction is the same with and without channeling at less than la JI^. 

recoiling ions, blocking effects (neglected in the DAMA calculation) are important. In fact in a 
perfect lattice no recoiling ion would be channeled (because of what Lindhard called the "rule of 
reversibility"), but due to lattice vibrations the collision with a WIMP may happen while an atom 
is displaced with respect to the string or plane where it belongs. If it is initially far enough it 
may be channeled. For an axial channel, the probability distribution function of the perpendic- 
ular distance to the row of the colliding atom due to thermal vibrations can be represented by a 
Gaussian, g{p) = {p/u\)e^^P~^^''^\ Thus, in our model the probability that an ion is channeled is 
given by the fraction of nuclei which can be found at a distance larger than a minimum distance 
P(,min (determined by p^. and the initial recoil angle) from the string at the moment of collision, 
Pch = Ip-^i^^^s{p) = e'^''""'"/^"'^. Notice that any uncertainty in is exponentially enhanced in 
the channeling fraction. Similar equations apply to planar channels. To obtain the total geometrical 
channeling fraction, we average the channeling probability over initial recoil directions, assuming 
they are isotropically distributed. We do it numerically by performing a Riemann sum once the 
sphere of directions has been divided using a Hierarchical Equal Area iso-Latitude Pixelization 
(HEALPix) [|T|] (a novel use for HEALPix). 

The channeling fractions we obtain are strongly T dependent (see e.g. Figs. 2,3 and 4.a). 
As T increases the probability of finding atoms far from their equilibrium lattice sites increases, 
which increases the channeling fractions, but the critical distances p^ become larger (~ cui at large 
enough energies) which decreases the channeling fractions. The fractions are smaller for larger 
values of c, as can be seen for Si and Ge in Fig. 3 and 4. We expect the fractions to be in between 
those for c = 1 and those for c = 2. Moreover, we have not included any dechanneling effects due 
to impurities or dopants. So the fractions we produced must be considered upper bounds. This 
is even more so for Nal and Csl, for which we have not found data equivalent to those provided 
by Hobler for Si (useful for Ge too, which has the same structure of Si) so we can provide up- 
per bounds on the channeling fractions instead of estimates, even forgetting about dechanneling 
effects. Upper bounds on channeling fractions in Nal are shown in Fig. 4.a for c = 1 and several 
temperatures. Using the fractions in Fig. 5. a corresponding to the temperature of the DAMA/Nal 
and the DAMA/LIBRA experiments, T= 293 K, the difference in the WIMP region corresponding 
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to the DAM A annual modulation signal differs only at 7 a is channeling is included or not included 
(as can be seen in Fig. 5.b taken from Ref. [12]). 

To conclude, the effect of blocking is important to understand the channeling of ions ejected 
from lattice sites by interactions with dark matter particles. Analytic models give good qualita- 
tive results but channeling data and or simulations are necessary to get reliable quantitative results. 
Finally, the daily modulation of the scintillation or ionization in crystalline detectors due to chan- 
neling, a potentially background free dark matter signal, deserves further investigation. 
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